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Abstract 

VLBI observations carried out by global networks provide the most accurate values 
of the precession-nutation angles determining the position of the celestial pole; as a 
rule, these results become available two to four weeks after the observations. Therefore, 
numerous applications, such as satellite navigation systems, operational determination 
of Universal Time, and space navigation, use predictions of the coordinates of the ce- 
lestial pole. In connection with this, the accuracy of predictions of the precession- 
nutation angles based on observational data obtained over the last three years is ana- 
lyzed for the first time, using three empiric nutation models — namely, those developed 
at the US Naval Observatory, the Paris Observatory, and the Pulkovo Observatory. 
This analysis shows that the last model has the best of accuracy in predicting the 
coordinates of the celestial pole. The rms error for a one-month prediction proposed 
by this model is below 100 microarcsecond. 

1 Introduction 

A prediction of the Earth orientation parameters (EOP) — the coordinates of the Earth's 
pole, Universal Time, and the coordinates of the celestial pole — is an extrapolation of EOP 
observations over a given time interval (the prediction length) following the last date of ob- 
servations. Such predictions are necessary for numerous scientific and practical applications, 
and numerous studies (see, for example, [1-6] and references therein) have been devoted to 
improvements in methods and analyses of accuracies of such predictions. Almost all studies 
in this field are concerned with predicting the coordinates of the Earth's pole and Universal 
Time. 

The reason for this is that, in contrast to these, precession-nutation models can be con- 
structed for the motion of the celestial pole (CP), with the accuracies of these models sat- 
isfying those required for numerous practical applications. In this situation, improvements 
in the accuracy of the precession-nutation theory simultaneously ensures improvements in 
the predicting accuracy. The last theory adopted by the International Astronomical Union 
(lAU) as the international standard, which incorporates the nutation theory IAU2000A (the 
lAU resolution B1.6 2000) and the precession theory P03 (the lAU resolution Bl 2006), has 
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an accuracy of 100 microarcsecond (/xas). However, this accuracy does not reflect tfie real 
error in the CP position calculated by the lAU theory, since the latter describes only part 
of the total motion of the Earth's rotational axis. In addition to the terms described by the 
theory, this motion also includes a number of terms which cannot be modeled at present, 
since they involve certain unpredictable effects, with the total contribution of these reaching 
400 yuas. The most important of these effects arc the free core nutation (FCN) and the trend. 

Thus, applications, which require only about 500 /ias accuracy in the precession- nut at ion 
modeling do not require any special predictions. However, applications with higher require- 
ments for modeling of the real-time CP motion, and all the moreso the CP motion prediction 
in the future, require modeling and predicting of the difference between the real CP motion 
and the lAU theory, called the celestial pole offset (CPO). The CPO values dX and dY 
are determined by Very Long Baseline Interferometry (VLBl) observations and the EOP 
services. There are several empiric models in use, which are continuously supported and can 
be used in practice for these purposes [7]. 

For each model, a CPO series is regularly pubhshed for each given time together with a 
prediction. However, only a few studies have considered the errors of the CPO predictions 
[8-10], and these studies have used historic data on the CP motion, which can distort the 
prediction accuracies, especially for short- term predictions [11]. In addition, only the study 
of Lambert [10] considers a model that is available for practical apphcations. Thus, users of 
predictions of the CP motion have no opportunity to obtain statistically based estimates of 
the associated errors. The only known attempt to estimate the real error of a CPO prediction 
was undertaken by the International Earth Rotation and Reference Frames Service (lERS) 
[12]. Unfortunately, this estimate analyzes only two predictions using unaccessible data and 
an unknown model. The graphical data presented in [12] indicate a CPO prediction accuracy 
of about 150 lias for prediction lengths of less than one year. 

Our present study analyzes real CPO predictions made by various EOP services in 2007- 
2009, using models saved in our database. These predictions are compared with the final 
CPO values determined by the International VLSI Service for Geodesy and Astrometry (IVS) 
using VLBI observations carried out by global networks [13, 14]. The statistical processing of 
these data provides rehable estimates for the quality of the CP prediction methods presently 
used. These estimates are calculated for the rms error, maximum error, and prediction 
efficiency. The first two quantities characterize the absolute error of the prediction, while 
the last estimates the relative error. 

2 Celestial pole offset determined from VLBI observa- 
tions 

It is well known that VLBI carried out on global networks provides the most accurate data 
on the CP motion. These IVS-coordinated observations require the participation of many 
observatories located in various countries [13]. In general, high-accuracy EOP data can 
be obtained from any set of VLBI stations with the necessary geometric parameters, most 
importantly, their size [15]. Since the data processing is very labor-intensive (we will analyze 
this below), VLBI observations become available with a delay reaching several months, which 
is insufficient for practical EOP applications, although the VLBI method is fundamental for 
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determining the CP position and Universal Time. To improve the immediacy of the EOP 
monitoring, IVS observations are incorporated in two special programs called Rl (as a rule, 
carried out from Monday 17^^ UTC to Tuesday 17^^ UTC) and R4 (as a rule, carried out from 
Thursday 18'*30"* UTC to Friday IS'^SO"' UTC). The observations incorporated in these 
programs have priority at all stages of the processing, reducing the delay in obtaining the 
EOP as much as possible. In addition, these observations are optimally scheduled to reach 
a higher EOP accuracy. 

The main principle of VLBI is that each station (radio telescope) of the VLBI network 
detects the radio-source signal independently. The resulting data are then sent to a corre- 
lator for processing. The correlation results include the interference delays and frequencies 
recorded in the form of IVS Mark-3 DBH (database v. 1 in the international terminology) 
standard files for the S and X bands. Further is the preliminary processing of the data, 
which takes into account the influence of the ionosphere, as well as data on the cable delay 
and meteorological and some other parameters both immediately detected at the stations 
and calculated theoretically. This results in final observational files (database v. 3 or 4 in 
the international terminology), which are used by the IVS centers for their calculations of 
the Earth rotation parameters. After processing the observations, the individual centers 
calculate IVS combined EOP series. All this data processing contributes to the time delay 
before the CPO data are generally available. Here, the time delay is the difference between 
the date when the EOP data become available and the epoch when the data are taken. 

Let us consider the approximately distribution of this time delay over the various pro- 
cessing steps for the EOP data determined in the IVS Rl and R4 programs. 

1. Data delivery to the correlator. Most VLBI observations carried out on global networks 
are recorded on magnetic disks and delivered to the correlator by express mail, which 
usually takes from three to seven days. This usually delays the start of the correlation 
processing by six to eight days. 

2. Data correlation. The correlation usually takes from two to five days. This also 
includes the preliminary analysis and data calibration, as well as the formation of the 
final observational file and its installation on the IVS servers for subsequent use at all 
processing centers. 

3. EOP calculations at the IVS processing centers. As a rule, first results become available 
at the centers forming the final observational files almost without any delay. At the 
other data-processing centers, the EOP data usually become available in one to two 
days^ . 

4. Calculation of the IVS combined EOP series. This usually takes from two to ten days. 

Our monitoring shows that the minimum delay for the IVS EOP series in 2009 is 12 
days and the maximum 37 days, with the median delay being 19 days. Thus, the CPO 
data of individual processing centers usually become openly available 1.5-2 weeks after the 
observations, while the IVS EOP series become available 2-4 weeks after the observations. 

^http://vlbi. geod.uni-bonn.de/IVS-AC/data/timeliness_2. html 
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The time delays for the observational data obtained in other IVS programs range from 
several weeks to several months, although these are almost insignificant for the final CPO 
data. Consequently, a CPO prediction for five weeks satisfies most practical applications. 
However, some other practical and research applications require accurate CPO predictions 
over longer time intervals. 

3 CPO models used 

To analyze the accuracy of the CP motion predictions, we used models freely accessible for 
research and practical applications. There are two types of such models, namely, the FCN 
models and the CPO models, which are related as follows: 

CPO = FCN + Y,Ti + Y,Pj, (1) 

where J2 Ti is the sum of trends, for example, resulting from errors in the precession model, 
and J2 Pj is the sum of (quasi) harmonic terms, for example, resulting from errors in the 
nutation model or geophysical processes. The FCN models are quasi-harmonic with variable 
amplitudes and phases [7], and their maximum amplitude reaches about 350 /xas. At present, 
the contribution of the second term reaches 200 //as and is continuously increasing; the 
contribution of the third term is probably below 100 //as. We have already noted that the 
sum of the CP coordinates determined by both the lAU precession-nutation theory and the 
celestial pole offset provides the most accurate data on the CP motion. Using the FCN 
models instead of the CPO models results in appreciable errors in the CP position, with 
these errors increasing with time. 

At present, a search for the most accurate CP coordinates and their predictions finds the 
following three CPO options. 

The NEOS CPO series. This series is calculated at the US Naval Observatory in the 
framework of the national US EOP service, and is an official combined online lERS 
EOP solution that includes, among the other data, the CPO series with a 90-day 
prediction. The algorithm used calculates the summary series of the CP coordinates 
using several selected VLBl EOP series and compares the results with the precession- 
nutation model [16, p. 73]. Formally, the CPO series is updated everyday; in fact, the 
new data become available when each new VLBI session is processed, i.e., two to three 
times per week, on average. 

The SL (S. Lambert) FCN model. This model is calculated at the Paris Observatory and 
is recommended by the lERS Conventions (2003) which collect the most up-to-date 
astronomical and geophysical models used for the reduction of astrometric observations 
[17, Ch. 5]. The combined lERS EOP scries calculated at the Paris Observatory and 
the analysis of its deviation from the precession-nutation theory is the basis for the 
FCN calculations. The FCN parameters are calculated via least-squares fits in shding 
two-year intervals, with the mean CPO values being removed within each interval 
[10]. In contrast to other series, this series contains no trends. The FCN parameters 
obtained within the last interval are used for a one-year prediction. Since 2009 the SL 
series is updated once per year — on July 1 (Lambert, private communication). 
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Table 1: Main CPO model characteristics. 



Model 


Type 


Forecast 


Smoothing 


URL 


USNO 


CPO 


3 months 


weak 


http:/ / www.usno.navy.mil/USNO / earth-orientation 


SL 


FCN 


1 year 


strong 


htt p : / /syrte . obspm . fr / ~ lambert / fen 


ZM2 


CPO 


2 ycrirs 


nioderMle 


hi Ip: / / www. gMo.spl).ru/cugiish/'Ms/persac 



The ZM2 CPO series proposed by the author. This CPO series is calculated at the Pulkovo 
Observatory using the IVS CPO series [7], and includes a two-year prediction. The 
construction of this series is described below. Starting from 2008, this series has 
been formally updated every day (it was updated about once per month before 2008), 
however, in practice, the new data become available only after the IVS series has been 
updated. 

The Table 1 presents the main characteristics of the given models. 

The ZM2 model is very simple and is calculated using Gaussian smoothing (high-frequency 
Gaussian filtering) of the IVS series individually for dX and dY: 

qi-eM-it-Ur/2a'], (2) 

Z^l PiQi 

where Xi are the original data containing n points detected with weights Pi at the time tj, 
x*{t) are the smoothed x values for the time t, and a is a smoothing parameter. When a 
increases, the smoothing also increases, i.e. the final series becomes smoother. Since the 
time t is arbitrary, this method can also be used for interpolation; see below. This feature of 
the Gaussian smoothing, as well as its insensitivity to coincident or very close ti, make this 
method preferable over the Whittaker and Vondrak methods frequently used in astrometry 
[18]. Autoregression fitting is used to provide the prediction for the ZM2 series. 

Figure 1 compares the three models with the IVS CPO series, and Fig. 2 presents the 
differences between these. The NEOS and ZM2 models are systematically close to each other, 
since these are constructed similarly, although the ZM2 data are smoother, as seems to be 
more natural. We suggest that the NEOS model is distorted by observational noise. Some 
systematic differences between these two models can probably be explained by differences in 
the VLBI data used, which have been mentioned above. We note also that the predictions 
obtained by these two models are very different (see the data shown in Fig. 1 after the last 
IVS point). 

In contrast to CPO models, the FCN SL model includes only the period of 430.2 days 
[10]. Therefore, in Fig. 1, the SL model is centered relative to the vertical axis, and shows a 
substantial difference between the model and observations (the IVS data). In addition, this 
model seems to be too smooth, and does not describe the Earth's rotation in detail. 

4 CPO prediction errors 

To estimate the accuracy of the CPO predictions, we used the predictions made at the 
US Naval Observatory (the NEOS model), Paris Observatory (the SL model), and Pulkovo 
Observatory (the ZM2 model) in 2007-2009, i.e., over the last three years. For the NEOS 
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Figure 1: The NEOS, SL, and ZM2 CPO models compared with the IVS data. 



and ZM2 models updated daily, we used a one-week interval for sampling the data. Thus, 
we used 153 NEOS, 121 ZM2, and 18 SL predictions. 

Every study of prediction errors includes an a posteriori analysis of the differences between 
the predictions and the final (observational) data. In our case, it is not entirely quite clear 
what data should be considered as final; both the CPO series obtained at the individual 
centers and the combined IVS or lERS EOF series can be used for this purpose. Wc decided 
to take the combined IVS EOF series as the most appropriate final data, since this series 
uses all the available VLBI data and avoids possible errors that arise in the construction of 
the combined lERS series. 

The large amount of noise in the observational data is the other difficulty in the analysis 
of prediction errors. This is true both for the individual and the combined CFO series, 
as we can see in Fig. 1. For the IVS series, the noise component of the CFO amplitude 
\/ dX'^ + estimated as the weighted Allan variation [19] reaches 118 /las for the total 
series obtained over 1984.0-2009.9 and 91 iias for the series over the last three years, i.e., for 
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Figure 2: Difference between the CPO models and the IVS data {dX is at the left, dY at 
the right; the NEOS, SL, and ZM2 models are shown in descending order). 

the data used here. 

Since the IVS data are given for the mean time of the daily observations, whereas all 
predictions are given for the beginning of the Greenwich day, we must reduce the IVS series 
to a standard smooth equidistant form with a step of one day, which we can do using 
an interpolation with smoothing, for example, Gaussian filter or spline smoothing. We 
used the first method applying formula (2). Here, we must choose the optimum smoothing 
parameter a. 

The difficulty is that the ZM2 and NEOS models use different, but quite similar, initial 
data and different smoothing parameters. Therefore, by varying the smoothing parameter 
for the IVS series when constructing the reference series for comparison, we can artificially 
bring this reference series closer to the series being compared, which will distort our estimates 
of the accuracy of the CPO predictions. 

To make our comparison as objective as possible, we calculated the prediction accuracy 
for six different smoothing parameters: a = 1, 2, 4, 8, 16, and 32. These calculations give 
similar results. Figure 3 presents the calculations for two smoothing parameters a=2, which 
is close to the NEOS series, and a=16, which is close to the ZM2 model. The original and 
smoothed IVS series are shown to the left, and the rms errors of the predictions to the 
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right. These errors are cahed external, since they are calculated relative to the "external" 
reference CPO series. The different smoothing parameters applied to the IVS series before 
comparison with the predictions result only in different smoothness of the curves, with no 
significant changes in the accuracy of the predictions obtained with the different methods. 
For completeness of the comparison. Fig. 4 shows the maximum prediction errors for a— 16. 
We arrive at the following conclusions about the accuracy of the CPO predictions: 

• the ZM2 model displays the highest prediction accuracy; 

• the NEOS model demonstrates a high accuracy for short-term forecasting, but the 
accuracy reduces with the prediction length; 

• the error of the SL CPO prediction is mainly determined by the systematic difference 
between the SL model and the IVS data. 

Note that, in contrast to the other two models, the Lambert model is an FCN model and 
not a CPO model, although it is recommended by the lERS Conventions as a CPO model 
(which is why we have included it in our comparison). For this reason, the SL model does 
not perform well compared to the NEOS and ZM2 models. Therefore, it is also of interest to 
study the accuracy of the predictions obtained using various methods for a given model. For 
this purpose, we repeated our analysis, but comparing the predictions with the final series 
obtained from the same model (instead of comparing with the final IVS series). The rms 
prediction errors for this case are shown in Fig. 5. These errors can be considered to be 
internal. 

To estimate the efficiency of the prediction methods, it was proposed in [1] to use the 
index P (predictability), defined as 

P(r) = l-^, (3) 

where r is the prediction length, crpir)^ the variance of the prediction errors for the given 
length, and the variance of the signal predicted. P — 1 indicates that the prediction 
coincides with the final results; this is typical for deterministic processes that are described 
completely by the model used. P < indicates that the prediction is completely inappropri- 
ate. In our case, the predicted signal is the CPO series. In general, the variance of the CPO 
series depends shghtly on the time interval considered. For the total IVS series obtained 
from 1984 to 2009, (Ts=0.18 milliarcsecond (mas) for dX and (75=0.23 mas for dY. For the 
three last years, as=0.20 mas for dX as=0.27 mas for dY. An increase in (Tg calculated 
over the last years is not surprising, since this variance is calculated for the original series 
including the trend that increases with time. For (7p(r), we use the rms prediction errors 
(see Fig. 3). Figure 6 presents the efficiency of the predictions. 

5 Summary 

We have analyzed the accuracy of CPO predictions calculated using three models that are 
currently available for practical use: the NEOS model calculated at the US Naval Obser- 
vatory, which is part of an online summary lERS EOP series; the SL model calculated at 
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Figure 3: External errors of the predictions obtained for various CPO models and two 
smoothing parameters for the IVS data. The original and smoothed IVS series are shown to 
the left, and the rms prediction errors to the right. 
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Figure 4: Maximum errors of the predictions obtained for various CPO models. 

the Paris Observatory and recommended by the lERS Conventions (2003) for high-accuracy 
modehng of the CP motion; and the ZM2 model calculated at the Pulkovo Observatory. 
The analysis used CPO predictions calculated in 2007-2009 and an a posteriori comparison 
of these predictions with the final IVS CPO series. Our statistical analysis of the differences 
obtained determined the following four characteristics of the accuracy of the prediction cal- 
culated using each method: 

• the rms error of the CPO prediction, 

• the maximum error of the CPO prediction, 

• the efficiency of the CPO prediction, 

• the rms error of the prediction compared with the same model. 

The first three errors describe the accuracy of the CPO prediction, while the last indicates 
the accuracy of the model self-reproduction. The ZM2 model demonstrates the highest CPO- 
prediction accuracy. In addition, only this model is efficient in forecasting both components 
of the CP motion, dX and dY ^ over a total prediction length reaching 100 days. The rms 
error of the one-month ZM2 prediction is about 75-100 yuas, depending on the smoothing 
parameter used. Since the celestial pole offset, which is the unpredicted part of the CP 
motion, reaches 400 //as and increases with time, this method improves the accuracy of the 
CPO prediction in ephemeris calculations by at least a factor of five. 

Only two of the models compared — the NEOS and ZM2 models — describe the CP motion 
completely, i.e. including the celestial pole offset. The model of Lambert is an FCN model; 
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Figure 5: Internal errors of the predictions obtained for various CPO models. 

i.e., it includes only one CPO component [see (1)], and cannot compete with the two other 
models in modeling the total CP motion, although it can be quite efficient in FCN predictions. 

Our main conclusions on the comparative accuracy of various predictions of the CP 
motion via a comparison of the rms prediction errors are supported by our analysis of the 
prediction efficiency (see Fig. 6). 

It is important that we have analyzed here predictions that are regularly calculated and 
freely accessible. Therefore, the results obtained can be helpful for both comparative studies 
of various predictions of the CP motion and estimates on the accuracy of the data obtained, 
as well as for choosing the most appropriate models for particular applications. 
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